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We examine the phase diagram of hadronic matter when the number of 
colours N c , as well as temperature and density, are varied. We show that 
in this regime a new percolation phase transition is possible, and examine 
the implications of this transition for extrapolations to physical QCD of 
the large- iV c limit. 

1. Introduction 

A powerful tool to investigate the features of non-perturbative QCD 
is the large- N c limit [I]: if the number of colours is sent to infinity in a 
controlled way, the theory undergoes huge simplifications and a number of 
results can be easily obtained. The underlying assumption that N c = 3 is 
"almost infinity" , however, seems to fail when it comes to predicting nuclear 
properties (as opposed, for instance, to hadronic ones). 

In the following we build a semiclassical model, based on percolation 
theory, to describe the possibility of a hidden phase transition along N c : 
its ingredients, the shape of a baryon and the likelihood for two quarks 
to perform some "exchange", can lead to a percolating system, meaning 
that, with a high enough density, we can have only local (i.e. confined) 
interactions yet achieve a global effect of correlations (exchange of quantum 
numbers) on a large scale. 

The final output of the model, a line separating ordinary confinement 
from a "confined percolating" phase, will then be compared to the decon- 
finement curves to assess the existence of a region, in the T-ijl-N c phase 
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space, where this regime is physically realised. A more thorough discussion 
on this model can be found in [2]. 



2. The percolation model 

We model dense baryonic matter by a regular cubic lattice with baryons 
sitting on the sites; each baryon will contain N c quarks, independently lo- 
cated according to a hard-sphere distribution with radius Aq<i D : 

/(x)«e(l-A QCD |x-x ccnter |) . (1) 

We parametrize density via the ratio e between the lattice spacing and 
twice the spheres' radius, so that, with ~p = Aq CD /8 (the reference case 
of [3] with baryons exactly touching each other), we now have p = Po e 3 - 
Thinking of the large- N c limit (where baryons are felt as quasi-static), we 
argue this description does not introduce large systematics at least down to 
N c ~ O(10), the coordination number of the 3D lattice considered [3]. 

The other important ingredient of the model is how quarks communicate 
with each other; this is encoded in a "squared propagator" function F(y), 
the probability of two quarks at distance y to exchange energy/momentum. 
As long as F scales as X/N c and has a sharp drop around a distance tt / Aq CD 
(tt ~ 1), we expect the result to be roughly independent of its exact shape; 
we consider three functions, corresponding to a step-function in coordinate 
space, a step-function in momentum space, and an exponentially screened 
version of the latterS 

Mv) = -tt ©f 1 - — t? — ) ; (2) 

Fk(v) = ^^tsm 2 ^— ) : Fs(y) = F K {y)e- M \y\ . (3) 
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With the baryon description and the "propagator" as input, together 
with the density e, we build the probability p that two baryons, with centres 
x^b, undergo some exchange on a short timescale ~ Aq,1 d as: 



p(N c 



J / A (x^)d 3 x A J / B (x B )d 3 x B (l-F(|x A -x B | 



(4) 



in the above formula (evaluated numerically for several choices of A, tt, N c G 
[2:80], e G [0.8:1.4] and the three F above), the quark-quark result is raised 



The maxima of Fk beyond the first one, albeit tiny, can trigger direct large-distance 
exchanges that we regard as nonphysical: in the following, we will rather consider Fs 
to get rid of this (high-density) effect. 
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to a power that can be either N c or N%, depending on whether we think of 
the exchanges as quark "flips" or "interactions" respectively. It has been 
seen in [3] that consistency with the large-iV c picture requires the latter 
choice, hence we set a = 2 from now on (a = 1 would give a decreasing 
p(N c ), see [2] for more discussion). We remark that the correlations among 
quarks, neglected in the above equation, were seen in [3] to amount to just 
a few percent correction. 

This work is an extension of [3j to a variable-density setting: in par- 
ticular, when e <C 1 baryons start to overlap, with the effect that limit- 
ing exchanges to nearest-neighbours would introduce systematic errors. For 
that reason, we implement a renormalisation-group-inspired approach, eval- 
uating a set of pi for nine relative distances (in lattice units) ranging from 
(1,0,0) to (2,2,2). These are then converted, with dedicated Monte Carlo 
simulations, into a probability Pb(Pi) for a cell of b 3 sites to achieve side-to- 
side percolation: the wrapping probability obtained in this way can finally 
be compared with the expected scaling, known from percolation theory [5j, 
to extract the threshold value ./V* at which percolation appears. As a con- 
sistency check, we verify substantial independence of the outcome of the cell 
size b in the range 3 < b < 7. 
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Fig. 1. Plots of N* as a function of the (adimensional) density, for a few sample 
choices of the parameters A, tt- Propagators are, from left to right: Ft, Fk, F$- 

The final result is a set of curves (for different propagators and different 
A,rr) in the form N*(jj) (Fig. [TJ: the translation of the density p into a 
physical baryonic density p raises, in principle, some ambiguities [2], but 
in practice the effect is negligible, so in this work they will be treated as 
coincident. 



3. Deconfinement line in the T-fi-N c space 

The next step is to locate the deconfinement line in the density- N c plane: 
above that line, indeed, baryons do not exist anymore, and the percolation 
model loses its validity. The deconfinement line, then, is essential in assess- 
ing whether the percolating phase is actually realised: in the following we 
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show how it is obtained, first at zero temperature and then at T > 0. 

The densities considered in the previous Section are the order of one 
baryon per baryonic size (i.e. ps ~ Aq CD /8): there, according to [6], mat- 
ter is still confined - at least for N c large enough. A parametric esti- 
mate of the zero-temperature deconfinement baryon chemical potential no 
can be obtained considering that the quark-hole screening (~ p^N c Nf, 
with Nf quark flavours) has to be the same order as the anti-screening 
gluon loop (~ N%): providing the correct dimensions with A QCD , we writqj 
Mo = N c ' Njr Aq CD ; neglecting three- (or more) body interactions, this 
iV c -scaling is quite a firm expectation [3J. On the other hand, the p = 
deconfinement point is iV c -independent, phenomenologically described by 
T c ~ (2/3)A QCD : we opt for a "minimal Ansatz" to parametrise the decon- 
finement line in the T — fx plane as 

2> 



fio J T c 2 A, 



QCD 

Employing the ideal gas formulas, we can express the relation (so far, at 
T ~ 0) between (baryonic) density ps and chemical potential /*b as 

^9f9s f°° j P 2 dp 

PB ~ 



(2nf 



I \ — n / ; P d = \i IP'B ->• -vb] \ (6) 

■A) [exp ±(^/ p 2 + m 2-p B ) +1 ' vj 'J 



(the second term accounting for antibaryons) , where gf is a flavour mul- 
tiplicity (one for Nf = 1, Nf(Nf — 1) for more), g s is a spin-multiplicity 
factor (g s = 2 for spin- 1/2 baryons), and m = N c Aq CD is the baryon mass. 
With a change of variables designed to put the iV c -scaling in evidence (note 
that we have automatically 7 = y/Nj), 



Mo 

we rewrite the density as 



^N c VN C A QCD N c 
7 = m ; a = V ; P = 7= > (?) 



_47rg I g 1 3 r f°° a 2 da "l / 8 n 



(2^)3^3/2^^^ 1 + exp[f3{ ^ T ^_ VW - c ^ )] 

This can be solved exactly at zero temperature (/3 — > 00), as shown in Fig. [2] 
for Nf = 1,3. Contrary to the findings in [3], we see that a percolating 
confined phase seems possible also in our three-colour world, but one has to 
push the density to extremely high values: this, however, has to be taken 
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Fig. 2. Log-log plot of the T = deconfmcment curves in the p-N c plane from 
Eq. [5] (and their large- N c limit) for Nf = 1,3, compared to the percolation curves 
obtained as described in Section [5] for the propagators Ft and Fs and some repre- 
sentative choices of the associated parameters. The horizontal line marks N c = 3. 
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Fig. 3. Density curves at deconfinement in the p-N c plane, for various temperatures 
T = 8T C (datapoints; the lines are only a guide to the eye); the decreasing lines 
with no datapoints are paramctrisations of two reference percolation curves for the 
two propagators Ft and Fs- to the pure mesonic contribution. Plots on the left 
(right) are for N f = 1 (3). 

with a grain of salt, the parametric assumptions made here being somewhat 
questionable towards N c = 3. 

In order to extend the calculation to nonzero temperature, we will pro- 
ceed numerically: here, there will be an antibaryonic contribution and, most 
important, we will not neglect the presence of higher-spin baryons. So, be- 
sides trading f3 for 6 using Eq. we consider that at a given N c = 2Q + 1 
we will have baryonic spins (rj + l)/2, with r] = 0, 1, . . . Q, associated to an 
additional energy r]A QCD /N c and with a spin-degeneracy of 2rj + 2: we then 
replace, in Eq. [SJ g s — > Y^rjfiv + 2) and add the cost of the spin-flip to the 
exponent in the denominator. The resulting curves in the ~p-N c plane are 



2 Unless specified by the subscript "q" , all chemical potentials are baryonic. 
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shown in Fig. O 

Looking at Figs. [2] and [3l it is clear that in the T-fx-N c space there 
is, besides the usual confined and deconfined phases, also a "confined but 
percolating" regime (upper right region in the plots). The question whether 
it actually extends down to N c = 3 is obviously delicate (the roughness 
of this model introduces systematics uncertainties, and it is not entirely 
clear if our world should be treated as a Nf = 3 or Nf = 2 setting, for 
instance), but if there is the possibility, it certainly has to be investigated 
experimentally at densities at least several times that of ordinary baryonic 
matter and at as low a temperature as possible. An important cross-check, 
at least ideally given the massive numerical effort involved, would be to 
perform lattice experiments at high N c . 

4. Conclusions 

The percolation model presented here, although rather qualitative, shows 
the possible presence of a phase characterised by confinement and, simul- 
taneously, large-scale correlations: in [3j it was argued that this could be 
identified with the "quarkyonic" phase proposed in [6]. In this work we ex- 
amine the shape, in the T-jj,-N c space, of the region where this regime occurs, 
finding that, provided the densities are high enough, there is a substantial 
possibility to observe its features in our world too. Based on the analogy 
between the percolating phase in hadronic matter and the metal-insulator 
transition, possible experimental signatures for this phase can be suggested 
[2]. A caveat is that, since the model is parametrically constructed in the 
large- N c spirit, its validity at N c = 3 could be hampered by deviations from 
the expected asymptotic scaling: at N c = 3 nuclear matter is a glass rather 
than a crystal, and nucleons are not fully classical objects. 
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